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' The conversion of neutron matter to strange matter in a neutron star 

, have been studied as a two step process. In the first step, the nuclear 

matter gets converted to two flavour quark matter. The conversion of 
two flavour to three flavour strange matter takes place in the second step. 
. The flrst process is analysed with the help of equations of state and hydro- 

CL(' dynamical equations, whereas, in the second process, non-leptonic weak 

JL I interaction plays the main role. Velocities and the time of travel through 

. the star of these two conversion fronts have been analysed and compared. 

I. Introduction 

^ , It has been conjectured that strange quark matter, consisting of almost equal 

' numbers of u, d and s quarks, may be the true ground state of strongly in- 

teracting matter |2] at high density and/or temperature. This conjecture 
is supported by bag model calculations for certain range of values for the 
strange quark mass and the strong coupling constant. By considering realistic 
values for the strange quark mass (150 - 200 MeV jll), it may be shown that 
the strangeness fraction in a chemically equilibrated quark matter is close to 
unity for large baryon densities. Such bulk quark matter would be referred to 
as "strange quark matter (SQM)" in what follows. 

The above hypothesis may lead to important consequences both for labo- 
ratory experiments as well as for astrophysical observations. Normal nuclear 
matter at high enough density and/or temperature, would be unstable against 
conversion to two flavour quark matter. The two flavour quark matter would 
be metastable and would eventually decay to SQM, releasing a finite amount 
of energy in the process. Such conversion may take place in the interior of a 
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neutron star where the densities can be as high as (8-10)po with po being the 
nuclear matter density at saturation OE]. If Witten's conjecture [T] is correct, 
the whole neutron star may convert to a strange star with a significant fraction 
of strange quarks in it. (Neutron star may also become a hybrid star with a 
core of SQM in case the entire star is not converted to a strange star. Such 
a hybrid star would have a mixed phase region consisting of both the quark 
matter and the hadronic matter [J].) Hadron to quark phase transition inside a 
compact star may also yield observable signatures in the form of Quasi-Periodic 
Oscillations (QPO) and the Gamma ray bursts (HUHl- 

There are several ways in which the conversion may be triggered at the 
centre of the star. A few possible mechanisms for the production of SQM in 
a neutron star have been discussed by Alcock et al ^OI- The conversion from 
hadron matter to quark matter is expected to start as the star comes in contact 
with a seed of external strange quark nugget. Such a seed would then grow by 
'eating up' baryons in the hadronic matter during its travel to the centre of the 
star, thus converting the neutron star either to a strange star or a hybrid star. 
Another mechanism for the initiation of the conversion process was given by 
Glendenning [3 . It was suggested there that a sudden spin down of the star 
may increase the density at its core, thereby triggering the conversion process 
spontaneously. 

Conversion of neutron matter to strange matter has been studied by sev- 
eral authors. Olinto ^I] viewed the conversion process to proceed via weak 
interactions as a propagating slow-combustion (i.e. a defiagration) front and 
derived the velocity of such a front. Olsen and Madsen and Heiselberg et 
al estimated the speed of such conversion front to range between 10 m/s to 
100 km/s. The combustive conversion front was assumed to have a microscopic 
width of a few tens to a few hundreds of fm in these calculations. 

Collins and Perry ^3], on the other hand, assumed that the hadronic matter 
gets converted first to a two fiavour quark matter that eventually decays to 
a three fiavour strange matter through weak interactions. Lugones et al 
argued that the hadron to SQM conversion process may rather proceed as a 
detonation than as a defiagration even in the case of strangeness production 
occuring through seeding mechanisms jlOj . 

Horvarth and Benvenuto examined the hydrodynamic stability of the 
combustive conversion in a non-relativistic framework. These authors inferred 
that a convective instability may increase the velocity of the defiagration front, 
so that a transition from slow combustion to detonation may occur. They 
argued that such a detonation may as well be responsible for the type II su- 
pernova explosions |1 7). In a relativistic framework, Cho et al [18\ examined the 
conservation of the energy-momentum and the baryonic density fiux across the 
conversion front. Using Bethe- Johnson JHj and Fermi-Dirac neutron gas ^Sl 
equations of state (EOS) for the nuclear matter (NM) and the Bag model for 
SQM, they found that the conversion process was never a detonation but a slow 
combustion only for some special cases. Recently, Tokareva et al modelled 
the hadron to SQM conversion process as a single step process. They argued 
that the mode of conversion would vary with temperature of SQM and with the 
value of bag constant in the Bag model EOS. Berezhiani et al Bombaci et 
al [2BI and Drago et al i24|, on the other hand, suggest that the formation of 
SQM may be delayed if the deconfinement process takes place through a first 
order transition so that the purely hadronic star can spend some time as a 
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metastable object. 

In this paper, we model the conversion of nuclear matter to SQM in a neutron 
star as occuring through a two step process. Deconfinement of nuclear matter 
to a two (up and down) flavour quark matter takes place in the first step in 
strong interaction time scale. The second step concerns with the generation of 
strange quarks from the excess of down quarks via a weak process. We may add 
here that this is the first instance where a realistic nuclear matter EOS is used 
to study the nuclear matter to SQM conversion as two step process. Drago et al 
[2i| . on the other hand, studied the burning of nuclear matter directly to SQM 
in detail by using the conservation conditions and the compact star models. 

To study the conversion of nuclear matter to a two-flavour quark matter, 
we here consider relativistic EOSs describing the forms of the matter in respec- 
tive phases. Along with such EOSs, we would also consider hydrodynamical 
equations depicting various conservation conditions to examine such conversion 
process in a compact neutron star. Development of the conversion front, as 
it propagates radially through the model star, would be examined. We would 
next study the conversion of two-flavour quark matter to a three-flavour SQM 
through a non-leptonic weak interaction process by assuming P equilibrium for 
the SQM. The paper is organised as follows. In section II, we discuss the EOSs 
used for the present work. In section III, we discuss the conversion to two- 
flavour quark matter. Conversion to three flavour SQM is discussed in section 
IV. In section V, we summarise the results. Conclusions that may be drawn 
from these results regarding the actual conversion process that may take place 
in a neutron star are also presented in this flnal section. 

II. The equation of state 

The nuclear matter EOS has been evaluated using the nonlinear Walecka model 
[2fi| . The Lagrangian density in this model is given by: 

£{x) = ^ 'iiji{i'j'^df_, -mi+ g^ia + .gt^itJ^7^ - gpip1^^^Ta)i!i - iw^^^w^^ 

i 

- ^brUnigaNO^)^ ~ ^C{gaN(^)'^ + 1pe{h''df_, - me)lpe (1) 

The Lagrangian in equation (1) includes nucleons (neutrons and protons), 
electrons, isoscalar scalar, isoscalar vector and isovector vector mesons denoted 
by ipi, ipe, c, uj^ and p"'''^, respectively. The Lagrangian also includes cubic and 
quartic self interaction terms of the a fleld. The parameters of the nonlinear 
Walecka model are meson-baryon coupling constants, meson masses and the 
coefficient of the cubic and quartic self interaction of the a mesons (b and c, 
respectively). The meson fleld interacts with the baryons through linear cou- 
pling. The u! and p meson masses have been chosen to be their physical masses. 
The rest of the parameters, namely, nucleon-meson coupling constant (^, ^ 
and ^) and the coefficient of cubic and quartic terms of the a meson self in- 
teraction (b and c, respectively) are determined by fitting the nuclear matter 
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saturation properties, namely, the binding energy /nucleon (-16 MeV), baryon 
density (po=0.17 fm~^), symmetry energy coefficient (32.5 MeV), Landau mass 
(0.83 TO„) and nuclear matter incompressibility (300 MeV). 

In the present paper, we first consider the conversion of nuclear matter, con- 
sisting of only nucleons {i.e. without hyperons) to a two fiavour quark matter. 
The final composition of the quark matter is determined from the nuclear mat- 
ter EOS by enforcing the baryon number conservation during the conversion 
process. That is, for every neutron two down and one up quarks and for every 
proton two up and one down quarks are produced, electron number being same 
in the two phases. While describing the state of matter for the quark phase 
we consider a range of values for the bag constant. Nuclear matter EOS is 
calculated at zero temperature, whereas, the two fiavour quark matter EOS is 
obtained both at zero temperature as well as at finite temperatures. 

III. Conversion to two flavour matter 

In this section we discuss the conversion of neutron proton (n-p) matter to two 
fiavour quark matter, consisting of u and d quarks along with electrons for the 
sake of ensuring charge neutrality. We heuristically assume the existence of a 
combustion front. Using the macroscopic conservation conditions, we examine 
the range of densities for which such a combustion front exists. We next study 
the outward propagation of this front through the model star by using the 
hydrodynamic {i.e. Euler) equation of motion and the equation of continuity 
for the energy density fiux ^jj- this study, we consider a non- rotating, 
spherically symmetric neutron star. The geometry of the problem effectively 
reduces to a one dimensional geometry for which radial distance from the centre 
of the model star is the only independent spatial variable of interest. 

Let us now consider the physical situation where a combustion front has 
been generated in the core of the Neutron star. This front propagates outwards 
through the neutron star with a certain hydrodynamic velocity, leaving behind 
a u-d-e matter. In the following, we denote all the physical quantities in the 
hadronic sector by subscript 1 and those in the quark sector by subscript 2. 

Condition for the existence of a combustion front is given by |28j 



where p is the pressure and X = {e + p)/nB'^, ub being the baryon density. 
Quantities on opposite sides of the front are related through the energy density, 
the momentum density and the baryon number density fiux conservation. In 
the rest frame of the combustion front, these conservation conditions can be 
written as pTHTzll^ : 



e2(p,X) < ei(p,X), 



(2) 



'^l^^lTl +Pl = ^2vl^2 + P2, 



(3) 



(4) 



and 



niViji — n2V2j2- 



(5) 
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In the above three conditions Vi (i=l,2) is the velocity, pi is the pressure, 

7i = . , is the Lorentz factor, oji — €{ + pi is the specific enthalpy and is 
V 

the energy density of the respective phases. 

Besides the conservation conditions given in Ij2l5|l , the condition of entropy 
increase across the front puts an additional constraint on the possibility of the 
existence of the combustion front. This entropy condition is given by [3(1] . 

siwi7i < S2W272 (6) 

with Si being the entropy density. 

The velocities of the matter in the two phases, given by equations ijSEJ, are 
written as |27] : 

„2 _ iP2 -Pi)(e2 +Pl) 



and 



(62 - ei)(ei +P2) 
{P2 -Pi)(ei +P2) 



(8) 



(£2 - ei)(e2 +Pi) ' 

It is possible to classify the various conversion mechanisms by comparing the 
velocities of the respective phases with the corresponding velocities of sound, 
denoted by Csi, in these phases. Thes conditions are 

strong detonation : vi > Csi, V2 < Cs2, 

Jouget detonation : vi > Cji, V2 — Cs2, 

supersonic or weak detonation: vi > Csi, V2 > Cs2, 

strong defiagration : vi < Csi, V2 > Cs2, 

Jouget deflagration : vi < Csi, W2 = Cs2, 

weak deflagration : vi < Csi, V2 < 0^2- 

For the conversion to be physically possible, velocities should satisfy an 
additional condition, namely, < < 1. We here flnd that the velocity 
condition, along with the eq.l(2j, puts severe constraint on the allowed equations 
of state. 

To examine the nature of the hydrodynamical front, arising from the neu- 
tron to two flavour quark matter conversion, we plot, in flg.l, the quantities 
vi,V2,Csi and Cs2 as functions of the baryon number density (ns). As men- 
tioned earlier, the u and d quark content in the quark phase is kept same as the 
one corresponding to the quark content of the nucleons in the hadronic phase. 
With these fixed densities of u and d quarks and electrons, the EOS of the 
two fiavour matter has been evaluated using the bag model prescription. We 
find that the energy condition (eqn.lQ) and velocity condition ( Vi'^ > ) both 
are satisfied only for a small window of « zLB.OMeV around the bag pressure 
51/4 = l60MeV. The constraint imposed by the above conditions results in 
the possibility of defiagration, detonation or supersonic front as shown in the 
figs.(l-3). 
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In fig.l, we considered both the phases to be at zero temperature. A possibil- 
ity, however, exist that a part of the internal energy is converted to heat energy, 
thereby increasing the temperature of the two flavour quark matter during the 
exothermic combustive conversion process. Instead of following the prescrip- 
tion for the estimation of temperature as given in references |El > study 
the changes in the properties of combustion with the temperature of the newly 
formed two flavour quark phase in the present paper. In flgures 2 and 3, we plot 
the variation of velocities with density at two different temperatures, namely, 
T = 50MeV and T = lOOMeV, respectively. These figures show that the range 
of values of baryon density, for which the fiow velocities are physical, increases 
with temperature. Figure 4 shows the variation of velocities with temperature 
for values of baryon number densities given by ub ~ 3po and 7po, respectively. 
In this figure, the difference between velocities vi and V2 increases with temper- 
ature of the two fiavour quark matter. In the present paper we have considered 
only the zero temperature nuclear matter EOS. On the other hand, equation 
of state of quark matter has a finite temperature dependence and hence the 
difference between vi and V2, varies with temperature. 

The preceding discussion is mainly a feasibility study for the possible gen- 
eration of the combustive phase transition front and its mode of propagation. 
Having explored such possibilities, we now study the evolution of the hydro- 
dynamical combustion front with position as well as time. This might give us 
some insight regarding the actual conversion of a hadronic star to a quark star 
and the time scale involved in such a process. To examine such an evolution, 
we move to a reference frame in which the nuclear matter is at rest. The speed 
of the combustion front in such a frame is given by Vf — —vi with vi being the 
velocity of the nuclear matter in the rest frame of the front. 

In the present work, we use the special relativistic formalism to study the 
evolution of combustion front as it moves outward in the radial direction inside 
the model neutron star. The relevant equations are the equation of continuity 
and the Euler's equation, that are given by |21j : 



and 



1 ,de de. 1 ,dv dv , ^ 
oj OT or dr Ot r 

where, w = ^ is the front velocity in the nuclear matter rest frame and k = ^ 
is taken as the square of the effective sound speed in the medium. 

Substituting these expressions for v and k in equations lO and itTIIll we get 

2v de 1 dv . 2^ 2ti „ , , 

OJ dr dr r 



n de ,^ 2 2v dv 



and 

uj dr dr 

Equations l(TT|l and l(T2| ultimately yield single differential equation that is 
written as: 

dv _ 2vkW^{l + v^) 
dr ~ r[4i;2 - fc(l + ^2)2] ' 
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The equation II 1311 is integrated, with respect to r(t), starting from the centre 
towards the surface of the star. The nuclear and quark matter EOS have been 
used to construct the static configuration of compact star, for different central 
densities, by using the standard Tolman-Oppenheimer-Volkoff (TOV) equations 
[32] . The velocity at the centre of the star should be zero from symmetry 
considerations. On the other hand, the 1/r dependence of the dv/dr, in ea. ljl3|l 
suggests a steep rise in velocity near the centre of the star. 

Our calculation proceeds as follows. We first construct the density profile 
of the star for a fixed central density. Equations Q and l|Hl then specify the 
respective fiow velocities vi and V2 of the nuclear and quark matter in the rest 
frame of the front, at a radius infinitesimally close to the centre of the star. 
This would give us the initial velocity of the front (— Wi), at that radius, in 
the nuclear matter rest frame. We next start with equation Ijl3ll from a point 
infinitesimally close to the centre of the star and integrate it outwards along 
the radius of the star. The solution gives us the variation of the velocity with 
the position as a function of time of arrival of the front, along the radius of the 
star. Using this velocity profile, we can calculate the time required to convert 
the whole star using the relation v = dr/dr. 

In figure 5, we show the variation of the velocity for values of the central 
baryon densities 3po, 4.5po and 7po, respectively. The respective initial velocities 
corresponding to such central densities are taken to be 0.43, 0.67 and 0.68. The 
figure shows that the velocity of the front, for all the central densities, shoots up 
near the centre and then saturates at a certain velocity for higher radius. Such 
a behaviour of velocity near the central point is apparent from the equation Ijl3|l 
above. The numerically obtained saturation velocity varies from 0.92 for central 
baryon density 3po to 0.98 for 7po- The existence of a saturation velocity, at large 
r, is apparent from the asymptotic behaviour of equation lfT3jl . A comparison 
with fig.l shows that for the densities 3po and 4.5po; the conversion starts as 
weak detonation and stays in the same mode throughout the star. On the 
other hand, for 7po, initial detonation front changes over to weak detonation 
and the velocity of front becomes almost 1 as it reaches the outer crust. The 
corresponding time taken by the combustion front to propagate inside the star 
is plotted against the radius in figure 6. The time taken by the front to travel 
the full length of the star is of the order of few milliseconds. According to 
the present model, the initial neutron star thus becomes a two fiavour quark 
star in about 10~^ sec. The results discussed above correspond to the case 
in which both nuclear as well as quark matter are at zero temperature. For 
finite temperature quark matter results vary only by a few percent of the front 
velocities for the quark matter at zero temperature. 

We would like to mention here that in the above discussions, the equations 
governing the conversion of nuclear to quark matter are purely hydrodynamic. 
There is no dissipative process, so that the combustion front continues to move 
with a finite velocity depending on the density profile. Furthermore, there is no 
reaction rate involved here as the deconfinement process occurs in the strong 
interaction time scale and hence can be taken to be instantaneous. This is 
certainly very different from the second step process, to be discussed in the next 
section, where the two fiavour matter converts to a three fiavour matter. Here, 
the governing rate equations are weak interaction rates which play a decisive 
role in the conversion. Comparing the total time ( = 10"'^ sec) taken by the 
combustion front to travel through the star with the weak interaction time scale 
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(10^^ — 10^^ sec), it is evident that the second step may start before the end 
of the first step process. In that case, perhaps, one should ideally considers two 
fronts, separated by a finite distance, moving inside the star i27|. In the present 
paper, we have taken a much simplified picture and considered the conversion of 
a chemically equihbrated two fiavour to three fiavour quark matter as the second 
step process. Our results may provide us with more information regarding the 
necessity of considering two fronts. 

IV. The conversion to three flavour SQM 

In this section we discuss the conversion of two fiavour quark matter to three 
fiavour SQM in a compact star. Similar to the discussion above, we assume the 
existence of a conversion front at the core of the star that propagates radially 
outward leaving behind the SQM as the combustion product. This conversion 
is governed by weak interactions that take place inside the star. 

For a three fiavour quark matter, the charge neutrality and the baryon num- 
ber conservation conditions yield 

Sub = Hu + Hd + Us (14) 

2n« = Ud + Us + 3ne (15) 

where rii is the number density of particle i (i= u, d, s and e). 
The weak reactions which govern the conversion of excess down quark to 
strange quark can be written as, 

d u + +Vi,~] s — > w + e~ + d + u^s + u (16) 

We assume that the neutrinos escape freely from the site of reaction and the 
temperature of the star remains constant. The non-leptonic weak interaction 
in such a case becomes the governing rate equation. The semi-leptonic weak 
decays, then, are solely responsible for the chemical equilibration which can be 
incorporated through the relations given below. 

Me- = - ^J.d = ^J's, (i7) 

where /i^ is the chemical potential of the i'th particle. The number den- 
sities {rii) of the quarks and electrons are related to their respective chemical 
potentials by. 



n,^9^ d^p/{2nf[f+ - n, (18) 
Jo 

where /"""and /~ are given by 



expHEp- ti)/T] + l ^ exp[{Ep + ^i)/T] + l■ ^^^^ 

In equations ifTsjl and ltT9|l . gi is the degeneracy factor and T the temperature. 
Eans. lll4ll9|l can be solved self consistently to calculate the number densities of 
quarks and electrons. 
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The conversion to SQM starts at the centre (r — 0) of the two flavour 
star. Assuming that the reaction region is much smaller than the size of the 
star, we have considered the front to be one dimensional. Moreover, as we are 
considering spherical static stars only, there is no angular dependence. The 
combustion front, therefore, moves radially towards the surface of the star. As 
the front moves outwards, excess d quarks get converted to s quark through 
the non leptonic weak process. The procedure employed in the present work is 
somehow similar to that of ref.^l], although the physical boundary conditions 
are different. 

We now define a quantity, 

air)^[nd{r)-ns{r)]/2nB (20) 

such that, a{r — 0) = qq at the core of the star. The quantity qq is the number 
density of the strange quarks, at the centre, for which the SQM is stable and 
its value Hes between and 1. For equal numbers of d and s quarks, a(r) — 0. 
Ideally at the centre of the star oq should be zero for strange quark mass nis = 0. 
Since s quark has a mass ~ ISOMeV, at the centre of the star oq would be 
a small finite number, depending on the EOS. The s quark density fraction, 
however, decreases along with the decrease of the baryon density towards the 
surface of the star, so that, a{r — > i?) — > 1 with R being the radius of the star. 
At any point along the radius, say r — ri, initial a(ri), before the arrival of the 
front, is decided by the initial two flavour quark matter EOS. The flnal a(ri), 
after the conversion is obtained from the equilibrium SQM EOS at the density 
corresponding to ri . 

The conversion to SQM occurs via decay of down quark to strange quark 
{u + d s + u) and the diffusion of the strange quark from across the front 
The corresponding rate of change of a{r) with time is governed by following two 
equations: 

f = -«(»)■ 

and 

§ - 

In the equation (I2TI1 R{a) is the rate of conversion of d to s quarks. Equation 
II22II yields the rate of change of a{r) due to diffusion of s quarks, with D being 
the diffusion constant. Following Olinto II], assuming the one dimensional 
steady state solution and using equations Ij21|l and[22lwe get: 

Da" - va - R{a) = 0, (23) 

where v is the velocity of the fluid. In equation lf23)l a — 

Conservation of baryon number flux at any position yields UqVq = n^Ws. The 
subscripts q and s denote the two flavour quark matter phase and SQM phase, 
respectively. The baryon flux conservation condition yields the initial boundary 
condition at any point r along the radius of the star: 

a'{r) = -^{a,ir)-af{r)), (24) 

where ai{r) and a/(r) are the values for the a{r) before and after the combustion, 
respectively. 
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The reaction rate for the non-leptonic weak interaction u + d ^ u + s is in 
general a five dimensional integral for non zero temperature and rris I^Mj . 
Here, instead, we have taken the zero temperature, small a limit lllj . 

i?(a) « ^GF^cos^9,s^n^^,^lu^^\ (25) 

lOTT 6 

where, Gp is the weak coupling constant and 9c is the Cabibbo angle. The 
above equation can be written in the following form: 

R{a) « — , (26) 

T 

where r = j^^^^^^G p'^ cos'^9cSin^9cfJ-u^ , here, depends on the position of the 
front. 

Following the line of arguments given in ref.[TT|, we write down the analytic 
expressions for D and v as: 

i?=.^^10-3(^)W/., (27) 
and 



-\/7 2(„.M-a,M) - <^«' 

Our calculation proceeds as follows. First we get the star characteristics for 
a fixed central baryon density pc- For a given pc, number densities of u, d and s 
quarks, in both the two and three fiavour sectors, are known at any point inside 
the star. That means ai(r) and a/(r) is fixed. Eans. (l25l28|l are then used to 
get the diffusion constant and hence the radial velocity of the front. 

The central baryon densities considered here are same as those of section III. 
Assuming that the neutrinos leave the star, the temperature is kept constant 
at some small temperature so that we can use the equation II25II . evaluated in 
the zero temperature limit. The variation of a(r) with the radius of the star is 
given in figure 7. The plot shows that a(r) increases radially outward, which 
corresponds to the fact that as density decreases radially, the number of excess 
down quark which is being converted to strange quark by weak interaction also 
decreases. Hence, it takes less time to reach a stable configuration and hence 
the front moves faster, as shown in the figs. 8 and 9. 

In Fig. 8, we have plotted the variation of velocity along the radius of the star. 
The velocity shows an increase as it reaches sufficiently low density and then 
drops to zero near the surface as d ^ s conversion rate becomes zero. Fig. 9 
shows the variation of time taken to reach a stable configuration at different 
radial position of the star. The total time needed for the conversion of the star, 
for different central densities, is of the order of 100 seconds, as can be seen from 
the figure. 



V. Summary and discussion 

We have studied the conversion of a neutron star to strange star. This conversion 
takes place in two stages. In the first stage a detonation wave is developed 
in the hadronic matter (containing neutrons, protons and electrons). We have 
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described this hadronic matter with a relativistic model. For such an equation of 
state the density profile of the star is obtained by solving Tolman-Openheimer- 
Volkoff equations. The corresponding quark matter equation of state is obtained 
by using the bag model. However, this quark matter equation of state is not 
equilibrated and contains two flavours. Matter velocities in the two media, as 
measured in the rest frame of the front, have been obtained using conservation 
conditions. These velocities have been compared with the sound velocity in 
both phases. 

For a particular density inside the star, flow velocities of the matter on the 
two sides of the front is now flxed. Starting from a point, inflnitesimally close to 
the centre, hydrodynamic equations are solved radially outward. The solution 
of the hydrodynamic equations gives the velocity profiles for different central 
densities. The velocity of the front shoots up very near to the core and then 
saturates at a value close to 1. The mode of combustion is found to be weak 
detonation for lower central densities. For higher central densities, the initial 
detonation becomes weak detonation as the the front moves radially outward 
inside the star. This result is different from that of ref.[21|, where the conversion 
process always correspond to a strong defiagration. The time required for the 
conversion of the neutron star to a two fiavour quark star is found to be of 
the order of few milliseconds. After this front passes through, leaving behind 
a two flavour matter a second front is generated. This second front converts 
the two flavour matter via weak interaction processes. The velocity of the front 
varies along the radius of the star. As the front moves out from the core to the 
crust, its velocity increases, implying faster conversion. The time for the second 
conversion to take place comes out to be ~ 100 seconds. This is comparable to 
the time scale obtained in ref.flT]. 

The comparison of the time of conversion from neutron star to two flavour 
quark star and the weak interaction time scale suggests that at some time during 
the passage of the first combustion front, the burning of two fiavour matter to 
strange matter should start. This means that at some point of time, there should 
be two fronts moving inside the star. On the other hand, our results show that, 
inside the model star, the burning of the nuclear matter to two fiavour quark 
matter takes much smaller time compared to the conversion from two fiavour 
quark matter to SQM. However, the consideration of two fronts might provide 
us with some more information regarding the conversion of neutron star to a 
final stable strange star. In the present case we have considered a two step 
process, there being only one type of front, inside the star, at any instant of 
time. 

Here we would also like to mention that ideally the second step should start 
with a non-equilibrated two fiavour matter [SSI ESI- Since this is a numerically 
involved calculation, in the present case we have taken the simplified picture of 
equilibrated quark matter. 

Finally, the burning of the nuclear matter to two fiavour quark matter is 
studies using special relativistic hydrodynamic equations. The actual calcula- 
tion should involve general relativity, taking into account the curvature of the 
front for the spherical star. We propose to explore all these detailed features in 
our subsequent papers. 
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Baiyon number density (fm-3) 



Figure 1: Variation of different velocities with baryon number density for T = 
MeV. The dark-shaded region correspond to deflagration, Hght-shaded region 
correspond to detonation and the unshaded region correspond to supersonic 
conversion processes. 




Baryon number density (fm-3) 



Figure 2: Variation of velocities with baryon number density for T = 50 MeV. 
Different regions correspond to different modes of conversion, where the nota- 
tions are the same as in fig 1. 
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Figure 3: Variation of velocities with baryon number density for T = 100 MeV. 
Different regions correspond to different modes of conversion, where the nota- 
tions are the same as in fig 1. 




Figure 4: Variation of velocities with temperature: (a)i;i for density 3po, (b) V2 
for density 3po, (c) Cgi, (d) Cs2 for 3po, (e) V2 for 7po, (f) ^^i for 7po and (g) Cs2 
for 7po. 
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Figure 5: Variation of velocity of the conversion front with radius of the star, 
for three different values of the central densities, namely, (a) 3po, (b) 4.5po and 
(c) 7po, respectively. Here po is the nuclear density. The initial velocity for the 
three cases are 0.66, 0.65 and 0.47, respectively. 



0.003 




8e+05 



Radius of the star (in cm) 

Figure 6: Variation of the time of arrival of the conversion front at a certain 
radial distance inside the star as a function of that radial distance from the 
centre of the star for three different central densities. Values of the central 
density corresponding to the curves (a), (b) and (c) are the same as shown in 
fig 5. 
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Figure 7: Variation of a(r) (as in the text) with radius of star, for different 
central densities, where, (a) corresponds to case for which the central density is 
3po, (b) for 4.5po and (c) for 7po- 




Figure 8: Variation of velocity of the two to three flavour quark conversion front 
with radius of the star for different central densities, where, (a) corresponds to 
case for which the central density is 3po, (b) for 4.5/9o and (c) for 7po- 
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radius of the star (in cm) 

Figure 9: Variation of time taken for the two to three flavour quark conversion 
front with radius of the star, for different central densities, where, (a) corre- 
sponds to case for which the central density is 3po, (b) for 4.5/9o and (c) for 

7 Po- 
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